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Comment on "Non-monotonicity in the 
Quantum- Classical Transition: Chaos Induced 
by Quantum Effects" 

In a recent Letter Kapulkin and Pattanayak pre- 
sented results regarding the chaotic behavior of the 
damped Duffing oscillator as it undergoes the transition 
from quantum to classical dynamics, this transition being 
induced by a model of continuous observation referred to 
as quantum-state diffusion They present evidence 
that a Duffing oscillator, sufficiently damped so that it is 
not classically chaotic, becomes chaotic in the transition 
region. If true, this would be a striking result. However, 
Kapulkin and Pattanayak did not calculate the Lyapunov 
exponent for the system, usually regarded as the litmus- 
test of chaos. Here we perform this calculation, which 
throws considerable doubt upon the conclusions in 

Since dynamical systems become very noisy as they 
pass through the transition, to calculate a Lyapunov ex- 
ponent at all, one must separate the sensitivity (chaos) 
induced by the deterministic dynamics from the unpre- 
dictability due to the noise. This can be achieved by com- 
paring the evolutions of the system for two nearby points 
in phase-space, using the same noise realization for both. 
If the separation of these trajectories, given by the usual 
distance in Hilbert space, is A(t), then the Lyapunov 
exponent is A = limt^oo {limA(o)^o (ln[A(t)/A(0)]/t)}. 

To calculate A numerically, A{t) must remain suffi- 
ciently small during the entire evolution that its dynam- 
ics is described to high accuracy by the dynamics of the 
system linearized about either of the trajectories. This 
can be realized by periodically re-scaling one of the tra- 
jectories towards the other along a line in the state-space 
of the system (which is in our case is simply the Hilbert 
space). This elegant procedure, a version of the Wolf 
method ^] , works because the linearized dynamics is un- 
affected by the size of A{t). We calculate the Lyapunov 
exponent by evolving a pair of trajectories for ~ 7000 
periods of the driving. To increase the accuracy we av- 
eraged the resulting exponent over 124 pairs, each with 
a different noise realization. For a given pair, the result 
obtained for the Lyapunov exponent, after evolving for 
a time t, initially falls as and then flattens out as it 
approaches and stabilizes at the true value. We report 
either this stabilized value, with the error determined by 
the spread over the different noise realizations (this error 
reduces as t increases), or, if the estimated value never 
flattens out, we report the resulting upper bound on the 
Lyapunov exponent. 

The Hamiltonian for the Duffing oscillator used in [l| 
is given by 
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FIG. 1: (Color online) The Lyapunov exponent for the 
damped quantum Duffing oscillator as a function of the in- 
verse system size, For = 0.01 the system is essentially 
classical. Red circles give the Lyapunov exponent, A, while 
blue triangles give upper bounds on A. (a) Small damping so 
that the oscillator is chaotic in the classical regime; (b) larger 
damping so that the classical oscillator has no chaos. 


tial factor of h does not affect the dynamics, the size 
(action) of the system with respect to h is entirely cap- 
tured by the parameter P = ^yh/{mu;P)^ where m and 
/ are, respectively, a mass and length scale. Damping 
given by P = —TP is induced by the quantum-state dif- 
fusion model [2I. Scahng time by Kapulkin and 
Pattanayak use the parameters Q = 1^ g = 0.3 and either 
r = 0.125 (classically chaotic) or F = 0.3 (classically 
unchaotic). Our results for these two cases are shown 
in Fig. 1. For the classically chaotic case the Lyapunov 
exponent decays as expected as h increases, and the sys- 
tem moves into the quantum regime. However, when 
the system is not classically chaotic, the Lyapunov expo- 
nent remains zero throughout. At least, our results place 
an upper bound on this exponent of 10~^. This throws 
considerable doubt on the conclusion in [1] that chaos 
emerges in the transition region for this system. 
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where X and P are position and momentum scaled so 
that [X, P] = and a; is a frequency. Since the ini- 
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